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ABSTRACT 

In this first year of our investigation we explored the role of compressibility and stratification 
in the dissipation of magnetic fields. The predictions of Mean Field Electrodynamics have been 
questioned because of the strong feedback of small scale magnetic structure on the velocity 
fields. In 2-D, this nonlinear feedback results in a lengthening of the turbulent decay time. In 
3-D alpha-quenching is predicted. Previous studies assumed a homogeneous fluid. This first year 
we present recent results from 2-D compressible MHD decay simulations in a highly stratified 
atmosphere that more closely resembles to solar convection zone. We have applied for NCCS 
T3E time to assist in the performance of our 3-D calculations. 


1. Introduction 

The goal of this project is to explore three key questions in solar dynamo theory: (1) In 2-D turbulent 
convection, how does the decay time of an imposed magnetic field vary as the field strength increases? (2) 
In a 3-D turbulent, compressible fluid, can a large scale, oscillatory field be maintained by a combination 
of small scale helical turbulence and a localized large scale shear, when feedback of the magnetic field 
on the momentum equation is included? (3) How much magnetic helicity is generated by a dynamo in a 
highly conducting fluid? How is the magnetic helicity distributed? When realistic boundary conditions are 
applied, what is the relative helicity in the ’’coronal field”? 

During the first year of this project we have made substantial progress on the first problem and have 
started work on the second. We have modified the existing compressible MHD code to include a random 
forcing by thermal blobs (see §2 below). This replaces the conventional forcing in phase space used in 
homogeneous turbulence models. The parameters of interest in these simulations are 7, the ratio of specific 
heats, A the dimensionless temperature difference across the layer, a the dimensionless viscosity, C the 
dimensionless resistivity, and Q the Chandrasekhar number. 


2. Equations and numerical methods 

Fully-compressible magnetohydrodynamics can be described by the nondimensionalized equations 

dp 

d t 


— u * V/> — pV ■ u, 


(i) 
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1 A-v M 2 

-u Vu- -Vp+ —Jjz + R^V-a- —*-{V 2 A)VA, (2) 

-u-VA + R~ 1 V 2 A ( 3 ) 

T(x,z,t) (4) 

together with boundary conditions at the sides, x = 0, A, and the top and bottom, z = 0,1, equation of 
state p = pT and prescription of the temperature T . Here the density p , velocity u, temperature T and 
magnetic flux function A (with the magnetic field B = — VA) have been nondimensionalized using the 
depth of the annulus d and the temperature T\ and density p\ at its top and the mean magnetic field B. 
Time is measured in units of the isothermal sound speed transit time at the top d/y/R+T\ where R+ is the 
gas constant. The resulting parameters 7, A, R e , Rm and M a denote, respectively, the ratio of the specific 
heats of the fluid, the dimensionless temperature difference across the height, the viscous and magnetic 
Reynolds numbers of the fluid and the dimensionless Alfven speed at the upper boundary. A denotes the 
horizontal aspect ratio of the computational domain relative to the height. The components of the viscous 
stress tensor a is defined by 

f $\i dxi 2 \ 

<?ij = + ~dx- ~ 1^'*^ u 7 with i = x ' z > 3 = x,z\ and l = x,z (5) 

In the absence of motion and current, equations (l)-(4) permit a static solution with 

T, = 1 + Az, p s =T T, P,=17 + \ A,=sm(*z)/7r (6) 


du 

~dt 

dA 

dt 


where m — 1/(7 — 1) • 


We assume that the motions are generated by randomly placed temperature sources which share a 
common amplitude and scale. These sources are defined using cubic B-splines where the temperature for a 
single source is 

ST{x, y,t;x 0 , z 0 , t 0 , t) = B s ^ B s (^~~j B » (“7“^) 

Here (i^,^,^) is the position and time of the temperature peak and \ an ^ r are the peak’s half width 
and lifetime. The continuous function B s (x) and its first and second derivatives vanish for |x| > 1 are 
continuous everywhere. It has the functional form of 


B s (z) 


2 (z + l) 3 , 

-6x 2 (x + 1) -h 1, 
< 6x 2 (x — 1) + 1, 
2(1 — x) 3 , 

l 0 , 


for — 1 < x < —0.5; 
for “0.5 < x < 0; 
for 0 < x < 0.5; 
for 0.5 < x < 1; 
for |x| > 1. 


( 8 ) 


The full temperature field has the form 


T(x,z,t) = T,(z)(\ +a'E?6T(x,y,t;x i ,z i ,ti,x,T )) 


( 9 ) 


where A^, a, x anc ^ T are specified for each simulation. The location of each fluctuation (xj,z*,ti) is chosen 
randomly in such a way that N sources are present at any time. Each source exists at a fixed location for 
its lifetime — r/2 <<<<,+ r/2. At the end of this interval, it is replaced by a new source at a new 
random position (x,*, z,*) with a new time center randomly chosen in the range t + r/2 < U < t + 3r/2 
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We assert that the viscous stress and vertical magnetic field vanish on the horizontal boundaries thus 
requiring the following for the vertical and radial velocity components and magnetic flux function: 

du d 4 

— = 0, w = gj = 0, at z = 0,1 (10) 

The vertical boundaries are taken to be periodic. In addition to these boundary conditions, we must satisfy 
the natural boundary condition on the pressure which arises through the combination of equation (2) and 
the boundary conditions (10). This condition states the the total normal force on the bounding surfaces is 
zero, and thus that the pressure gradient there must balance all other normal forces. 

The system of equations (1) through (4) are solved using a sixth-order compact difference scheme with 
fourth-order time accuracy (see Hurlburt 1997 J. Set. Comp, for details). For a 192 x 192 computational 
grid, each solution requires approximately ten seconds per step on an 180Mhz SGI R5000 Indy, with about 
270 steps required to advance the solution one time unit. 


3. Solutions 

We have begun by establishing the decay rate for essentially kinematic fields. Because we have 
introduced a non-standard forcing, it is important to understand the reference fields before proceeding with 
the parameter studies. Figures 1 & 2 show two of the reference solutions: figure 1 is a weak field case 
(plasma /3 = 1000) and figure 2 is a stronger field case (plasma f3 ~ 100). For both solutions R e — 12, 500, 
Rm = 50, 000, AT — 0.7, 7 = 1.4, and 9 = 0.5. Time traces of the magnetic energy are shown in figures 3. 
The straight solid line shows the Ohmic decay rate for the actual diffusivity. 

Figure 3 shows an example of the decay of the magnetic energy in a stratified layer for four values of 

for R e = 11, Rm — 22, 7 = 1.6666, A = 5, a = 0.7 and x — T “ 4. Kinematic solutions (dot-dash) with 
M a = 1.44 x 10 -3 are amplified by a factor of 10 and then decay rapidly. As the field strength increases 
over the range M a = 0.10, 0.18 and 0.32 (dash, dot and solid respectively), the field strength is amplified 
less and decay at the slower ohmic rate for longer intervals. The straight solid curve represents the ohmic 
diffusion rate. 

We presented preliminary solutions at the the Solar Physics Division Meeting in Montana. A paper on 
the 2-D results is in progress. 




